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Abstract 

Thermal properties of few-fermion (n < 5) systems are investigated. 
The dependence of the heat capacity on the topology and shape of the 
cavity containing the particles is analyzed. It is found that the maximum 
of the heat capacity, occuring at low T, discussed recently by Toutounji 
for a system with n = 1 fermions, is even more visible for n = 2, but fades 
away for n = 3 and 4. For large T, the classical behavior is obtained; 
however, when T — I 0, the heat capacity tends to zero exponentially, not 
linearly, as in macroscopic and even mesoscopic systems. The physical 
relevance of these results is discussed. 

1 Introduction 

Recent progress in nanophysics has stimulated the interest in studying particles 
confined to small volumes or in atypical geometries. If we are focused on the 
theoretical study of thermal properties of such systems, the contribution of 
translational degrees of freedom of the constituents to the partition function of 
the system might be evaluated taking into account the discrete character of the 
particle wave vectors. 

To be more specific, let us mention that, for non-interacting many-particle 
systems, the translational part of the canonical partition function is evaluated 
in the thermodynamic limit, as an integral over the particle momenta; however, 
for few-particle systems, confined to small cavities, this quantity should be ob- 
tained as a discrete sum over the quantized values of the particle wave vectors. 
Recently, Toutounji 1 analyzed the thermal properties of an electron confined 
to a segment, using the canonical ensemble, showing that the heat capacity of 
this one-particle system has a bump, at "low" temperatures, i.e. when fc^T is 
much lower than the ground state energy (the quotation marks intend to remind 
that, in this context, the room temperature is a "low" temperature). Of course, 
as T increases, the heat capacity reaches asymptotically the value fcs/2. The 
existence of this bump is smeared out in the thermodynamic limit. We can 
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presume that other interesting effects, specific to few-particle systems, are de- 
tectable only taking into account the discrete character of physical parameters 
describing such systems. A critical discussion of Toutounji's approach, from the 
perspective of the foundations of statistical mechanics, has been done by Lungu 

It is legitimate to wonder about the physical relevance of few-particle ther- 
modynamics. There are several aspects to be discussed: (1) is it meaningful to 
speak about thermal properties of few particle systems, knowing that thermal 
physics works essentially with large systems? (2) if the response is positive, is 
the canonical ensemble appropriate for such investigations? (3) there is a real 
physical interest for the study of thermal behavior of few-particle systems? 

If the particles are confined to a cavity in thermal equilibrium with a thermo- 
stat, they reach thermal equilibrium not through mutual interaction (the typical 
case studied by statistical mechanics) , but mainly through interaction with cav- 
ity's walls and with thermal radiation. In fact, the calculation of the one-particle 
partition sum for particles in simple external potentials is a common exercise 
in several books of quantum mechanics or quantum statistical physics (see for 
instance [3], [4]). However, as the thermodynamic relations are obtained, in 
general, in the thermodynamic limit, the significance of results referring to few 
particle systems, must be considered with caution [2J. 

Concerning the second issue, let us recall - invoking again the detailed analy- 
ses done by Lungu [2J - that the grand-canonical system is the most appropriate 
for statistical mechanics calculations, because the quantum statistical proper- 
ties of particles are automatically taken into account. However, we can use the 
canonical ensemble, if we correctly define the quantum states and their degen- 
eracies. This can be done quite easily for few particles, but becomes a quite 
cumbersome exercise if the number of particles increases. It is important to 
stress the fact that the canonical partition function of a quantum few-particle 
system does not factorize, in terms of one-particle partition functions. This is 
why the results obtained for a one-particle system have no relevance for several- 
particle systems. 

Finally, the study of few-particle systems is not an academic exercise, as such 
systems are available and have very interesting properties, the most popular 
being the electrons in quantum dots and atoms in Bose-Einstein condensates 
(BEC), see for instance [5], [6]). Investigations in quantum field theory provide 
another reason for calculating the one-particle partition function, for instance 
of a particle in an anharmonic potential [7]. Let us also mention that the last 
paper of Feynman presents a method for approximating the quantum partition 
function with a classical one [8] . 

The starting point of the investigations presented in this paper is the attempt 
to reconcile the heat capacity behavior obtained by Toutounji for a one-electron 
system - which shows an maximum at low temperature, as already mentioned 
- and the well-known aspect of this quantity, for a macroscopic body, which 
starts from zero at T = 0, increases smoothly while T increases, reaching an 
horizontal asymptote for T — > oo. In fact, the same smooth behavior occurs 
for mesoscopic systems with as few as n — 14 particles, as shown by Lungu 
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[2 . So, we are trying to bridge the results obtained by Toutounji, for a system 
with n = 1 particles, with the results obtained by Lungu, for a system with 
n = 14 particles. We evaluated approximately the canonic partition function, 
for a system of 2, 3 and 4 free fermions, taking into account the exclusion 
principle - so considering that the particles are indeed fermions - and found 
out that the low-T maximum occurring for an one-electron system decreases for 
systems with 3 and 4 particles (even if it is enhanced for a 2-particle system), 
and the aspect of the heat capacity plot becomes smoother. So, our results 
bridge the one-particle behavior of the heat capacity, obtained by Toutounji, 
with the mesoscopic behavior, obtained by Lungu, which in fact does not differ 
qualitatively from the macroscopic behavior. 

However, there is another "anomaly" which is noticed in the heat capacity of 
few-body systems. If, for macroscopic systems (see for instance [9], Ch.6), and 
even for mesoscopic ones [2J, the heat capacity decreases like T when T — > 0, in 
the cases discussed here (1, 2, 3 and 4 particles) it behaves like e~ 1 / T /r 2 , where r 
is a scaled temperature. Such a behavior is compatible with Toutounji's results, 
even if this issue is not explicitly discussed in his paper (see [T], Fig. 4). 

We also take advantage of the fact that the same - very simple - theoretical 
methods, used in Toutounji's paper, can be applied in order to point out the 
dependence of the heat capacity on the topology and shape of the system. After 
discussing the thermodynamic properties of a particle confined to a segment, we 
examine the similar problem, for the same particle confined to a circle obtained 
by bending that segment - so, an ID body having the same length, but a 
different topology. As the one-particle energy spectrum are different in the two 
cases, and as the heat capacity is sensitive to the aspect of the energy spectrum, 
we find that the heat capacity of a particle confined to a segment, and on a circle 
having the same length, are different. In this way, we obtain a simple example 
of topology-dependence of the heat capacity. Then, the same investigations are 
repeated for systems of 2, 3 and 4 particles, having the same confinement as 
mentioned before. 

Such a change of topology can be easily done, for charged particles moving in 
a thin toroidal cavity of length L: a bias applied on a small portion 6 of this torus 
will oblige the particles to move on the rest of the cavity, topologically equivalent 
(if the thickness of the cavity tends to zero) to a segment of length L — S. With 
a zero cost in energy, we changed the topology of the body and, implicitly, its 
thermal behavior. Another example is provided by the benzene molecule: if it 
is broken, the energy levels are quantized due to the impossibility of electrons 
to leave the linear molecule; by contrary, for electrons on the benzene ring, the 
quantization results from a cyclic condition [10] . The contribution of ir electrons 
to the heat capacity is different for a gas of intact benzene molecules, and for 
one of broken molecules (the ring can be broken, for instance, by ultraviolet 
irradiation). 

Another situation under scrutiny refers to few-particle systems confined to 
3D cavities: a rectangular (parallelepipedic) , cylindrical and spherical cavity. 
In the first case, the shape of the cavity can be easily modified, at constant 
volume, and the dependence of the heat capacity on the shape of the cavity can 
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be easily seen. 

The structure of this paper is the following. In Section [2] the partition 
sum, internal energy and heat capacity of a few-particle system, confined to a 
segment or to a circle of equal length, are obtained and compared. The values 
of the heat capacity differ drastically, in the two cases, for low temperatures, 
but tend to the same classical value, for high-T. In order to avoid irrelevant 
complications, scaled temperatures have been used. All the results are exact, 
being expressed in terms of the elliptic theta 3 Jacobi function #3 and of its 
first and second derivatives. The limiting cases of low- and high-T are obtained 
explicitly. Section [3] is devoted to a similar exercise, developed for particles 
confined to 3D cavities, as described in the previous paragraph. In the case 
of the spherical and cylindrical cavities, the results are no more exact, as the 
energy eigenvalues are obtained using an asymptotic approximation of the roots 
of Bessel functions; however, this approximation is quite accurate, even for the 
smallest values of the quantum number. Also, the summation involved in the 
evaluation of the partition sums must be made over two quantum numbers, 
running over N, resulting expressions similar to the Jacobi theta functions, but 
somewhat more complicated. However, the physical interpretation of results 
can be done without difficulty. In Section 31 we present final considerations and 
conclusions. 



2 A few-fermion system on a segment and on a 
circle 

A particle confined to a segment of length L is described by the Schrodinger 
equation for the infinite square well potential: 



V(x) = {^ X< \ (1) 
00 , x > 



2 

The impenetrability of walls means that: 



which constraints the wave vector to the values: 



nn , . 

k ns = — , n = ±l,±2,... (3) 

and consequently produce the energy eigenvalues: 

fc2p t2 2 

E ns = y^ = ±L_ n 2= Elsn i (4 ) 

2 m 2 m L 

where 
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is the ground state energy of the particle on the segment. 

The only effect of including the negative values of n is a phase factor multi- 
plying the wave function [TT], so they can be disregarded. 

If the particle moves on a circle of length L, the equation ^ is replaced by 
a cyclicity condition: 

ip (x) = i/j (x + L) (6) 

giving: 

Here and hereafter, the index s will be used for "segment", and c - for "circle". 
The energy levels are: 

h 2 k 2 2h 2 ir 2 



nc " 2 tp 2 



2m m L 2 

with 



n 



Ei c n' (8) 



E n+ i, c -E nc = _(2n + l), E n+1 , s - E ns = — (2n+l) (11) 



2ft 2 ^ 2 

Blc - 9 

m 

is the ground state energy of the particle on the circle. Let us notice that, as 
an effect of replacing the rigid- well condition @ with the cyclic condition ([5]), 
we have the relation: 

E lc = AE ls (10) 

As the energy of excited states is obtained from the ground state energy by the 
same relation, for the segment and for the circle (compare Eqs. ((4]), (Tl0| ). the 
energy levels on the circle are 4 times larger that the corresponding ones on the 
segment; this fact will have important consequences for the thermal behavior. 
The distance between consecutive levels: 

2h 2 v 2 

2m -L 

has of course the same property, reflected directly in the density of states and 
heat capacity. In fact, the most popular formula expressing the heat capacity 
in terms of the density of states g (E) is: 

C =Jf j ' Ef FTi (E,T)g(E) dE (12) 

It is properly invoked when the energy level are extremely dense, but it is evident 
that it shows that, at least at low T, a larger level spacing means a smaller heat 
capacity. 

As the energy levels are proportional to n 2 , the partition sum for the particle 
on a segment, Z s , and for an identical particle on a circle of the same length, 
Z c , can be expressed as: 

oo oo 

Z s = exp (-f3E ns ) = J2 if = Z (q s ) , (13a) 



ex P ( j2 ) =exp(-/3£ s i) ; (13b) 
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Z C = J2 ex P (-P E nc) = € =Z (q c ) , (14a) 

n=l 

g c = exp(-2/3— - — J =exp(-PEci) ■ (14b) 



We shall adopt a more general notation: 



where we have introduced the scaled temperature r: 

The sums in Eqs. (fl3|) - (|14p can be expressed in terms of the Jacobi theta 3 
function: 

oo 

tf 3 (0,<z) = l + 2^ fe2 = l + 2(q + q i + q 9 + ...) (17) 

fc=l 

So, the partition function is: 

Z(q) = ^ld 3 (0,q)-l} = q + q 4 + q 9 + ... (18) 

where g may be either q a or g c . 

Due to Eq. (fTDj) . we have the relation: 

9c = cxp (- /3 Eke) = exp (- 4 /3 E ls ) = q 4 s . (19) 

The thermodynamics can be obtained from Eqs. (|13|) - (|14p. using the well- 
known relations: 

<>-- B -w-"* a -w- « 

As usual, the cases of low and high temperatures are particularly important; in 
terms of g, the low-T limit corresponds to: 

g = expf-^ ~expf-^^ ~exp(-oo) = (22) 

and the sum l|17p is rapidly convergent. However, for high-T, 

9 = exp(-i)~exp(-i)=e° = l (23) 



G 



and the sum (jTTJ) converges very slowly. In this case, it is convenient to use the 
following relation, obtained from Wolfram - EllipticTheta3(09.03.06.0032.01) 



0s(O,g) =03 



^0,cxp (-^j) = Vkt$3 (0,exp(-7r 2 r)) , r>l. (24) 



It corresponds to eq. (8) of [T]. 

This means that, for high T, using Eqs. (fT8|) and ([24]) . the partition sum 
might be more conveniently written as: 



Z il) = \ [*» (0» ?) - !] = ^ ^3 (0, exp (-tt 2 t)) - l] , q = exp 

Of course, both Eqs. ([T£f and ([2"5j) are exact. 

In the limit of low T, and implicitly of small q: 



1 

(25) 



InZ (g) = lng + 9 3 - \ f + q* + ± g 9 - g 11 - ~ q 12 + q 14 



(26) 



and from the second equation in (|20[) we get: 

U = E x [1 + 3 q 3 - 3 <? 6 + 8 q H 

and 



t7 
Ex' 



-> 1 if r ->• 



(27) 



(28) 



As expected, at T = 0, the internal energy coincides with the ground state 
energy. Also, using the first equation in (|21[) : 



C = 9 fc B \ exp ( -- 



1 - 2 exp - 



64 

"77" exp 

9 \ T 



(29) 



so 



C -> if t — ^ , 



(30) 



as requested by the third principle of thermodynamics. In the r — > limit, the 
ratio C s /C c is: 

C.s 1 9^i s 



C c ^16 eXP A: B T 
At high T, again with the second equation of ((21 



(31) 



U = ^k B T + (rexp (-tt 2 t)) 



(32) 



and the heat capacity: 



(33) 
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so the Dulong-Petit law is obtained. As a trivial consequence, 



a 



-4- = l + 0(T 2 exp(-^ 2 r)) . 



(34) 



Comparing (pTTj) with (|34p . it is clear that, at low T, the quantum specificity of 
each system is exponentially dominant, and at high T - exponentially insignifi- 
cant. 

After getting some physical insight on the physical behavior of the systems, it 
is useful to obtain exact expressions for U and C . Let us introduce the following 
notations: 



,/,.,(!. , , = d ^j°' g) = 2 f; fc 2 q fc2 - 1 = 2 (1 + 4 q 3 + 9 g 8 



^'(0,?) = 



dg 

d 2 3 (0,g) 
dg 2 



fc=i 



= 2 ^ fc 2 (fc 2 - 1) q k ' 2 - 2 = 2 (12 q 2 + 72 g 7 



(35) 
(36) 



k=2 



warning about the fact that, traditionally, these symbols are attributed to the 
derivatives of $3 with respect to its first argument. 



U 



.E 



'TT T 

k 



U = E iq y : f , lowT,g~0 
i? 3 (0,g) - 1 

^ 1 f , c-- V ) - ttV 1 / 2 c-^ T d' 3 (o, e— V ) 



7TT z? 3 (0,e- 7r2 " r ) - 1 



(37) 



, high T, q ~ 1 



. (38) 

Both expressions for [7 are of course exact. Similarly, for the heat capacity, at 
low T it is convenient to use: 



C(q) =k B (lng) 2 



and at high T: 



■d" (q) 

[$3 (<z) — I] 2 ' ' fa (?)"!] '" ^ fa (?)-!] 



(39) 



d 1 r 3/2 ^ 3 _ 7r 2 r 5/2 £ -^t tf/ ^ g-^r^J 

c = V "d7 



/ttt i? 3 (Q,e-^ T ) - 1 



, high T, 9 ~ 1 . 

(40) 

Again, both Eqs. (|39|) and (l40l) are exact, but the series expansion in q of 
Eq. (|34j) . and the series expansion in r in EWq. (|40j) . are rapidly convergent, at 
low and high T, respectively. 

As we could see, the quantities T, m, L enter through the expression: 
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h 2 ir 2 
m L 2 



(3 



1 



m L 2 m i 2 T 



(41) 
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so the values of m, L separately are not relevant. Replacing the constants with 
their values and using atomic units: 



2^2 



2/3 



958. 34 



m L 2 Tm a L 2 ' 



(42) 



where m a , L a are expressed in amu and Angstrom. So, the only physical variable 
entering in the expression (|2"Tj) of the heat capacity is the scaled temperature r: 



q s = exp ^- /3 
q c = exp I -2/3 



ft 2 ^ 2 
2mL 2 

h 2 ir 2 
mL 2 



— exp 



exp 



Tm a L 2 a 
240 



Tm a L 2 a 
958.34 



(43) 



(44) 



Interesting effects are expected to occur at t < 1. Assuming that L a ~ 10 A, this 
means a temperature T < 2/m a ~ 1000 K for electrons and about 0.1 . . . 0.01 K 
for atoms. 

As the quantity q takes different values if the particle is confined to a segment 
or on a circle, according to Eq. (|19p. the heat capacity in these two cases is 
different. It depends on the topology of the body. 

Let us now investigate the situation of a system with more than one fermion, 
confined to a segment; in the case of the confinement on a circle, the results 
are absolutely similar. The energy levels of the system is obtaining filling up 
the one-particle levels and taking into account the degeneracy of states. For 
n = 2, 3, 4 particles, the partition functions are: 



Z 2 = exp 



3=1 

oo oo 



2 J1 

T 



Z 3 = 2j2J2 ex P 

j=i fc=i 



oo oo 

+ 4 EE ex p 

3=1 fc=l 

2f + (j+kf 



j 2 + (j + kf 



2 EE ex p 

3=1 fe=l 



(45) 



j 2 + 2(i + kf 



OO OO OO 



8 EEE-p 

j=l fe=l n =l 

2j 2 



j 2 + (j + k) 2 + {] + k + n) 2 



z * = E E ex p 

3=1 fe=l 



2(j + k) 



2 1 



OO OO OO 



OO oo oo 

4 EEE ex p 

oo oo oo 

4 EEE ex p 

j=l k=l p=l 

oo oo oo oo 

16 EEEE-p 

j—1 k—1 p—1 q—1 



4 EEE ex p 

3=1 fe=lp=l 

j 2 + 2(j + k) 2 + (j + k + p) 2 



(46) 



2j 2 + {] + k) 2 + {j + k+ P ) 2 



J 2 + (j + k) 2 + 2(j + k + p) 2 



(47) 



j 2 + (j + k) 2 + (j + k + p) 2 + {j + k+p + q) 



2 1 
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Clearly, each such expression has the form: 

(~ £ -f) +.-], (48) 



Z„ = cxp 



Sir, 

gon + gin exp ( — ] + g 2n exp 



where £ n is the Fermi energy of the n-particle system, e\ n , E2n, ■ ■ ■ ~ the excited 
states and gj n - their degeneracies of states. They are generalizations of the 
one-particle partition function: 

oo / -2 \ 

Z 1 = ]Texp i- 3 - . (49) 
j=i \ / 

As Z n for n > 1 cannot be expressed in terms of a finite number of known 
functions, for the evaluation of the heat capacity, the series in Eqs. (14"5"|) - (|47p 
will be cut at some value of I, fc, etc. So, it is useful to have an idea of errors 
introduced in this way. 

Let Z n (max = J) denote the finite sum obtained from Z n , if we keep only 
the terms of the sums with summation indices not larger than J. The only 
case when the finite sum can be compared with the exact value is n = 1. We 
find that InZi (max = 100) approximates \nZi with an error less than 10~ 15 
for r < 10 and In Z\ (max = 10) approximates In Z\ with an error less than 
10 -13 for r < 4 and less than 10~ 7 for r < 8. So, we can presume that, for 
a qualitative description of the heat capacity, it will be acceptable to use the 
sums Z n (max = 10). However, for n — 1, we have used Z n (max = 100); we 
avoided to work with the exact partition sum, as the exact calculation of the 
heat capacity involves the derivative of $3 with respect with its second argument, 
which is a very complicated function (not to be confused with EllipticTheta 
Prime 3) and is not very appropriate to be used in Mathematica calculations. 

It is convenient to use the scaled temperature r instead of the physical one 
T, and to define an effective heat capacity per particle: 



1 dU n 1 d ( 2 dlnZ„ 



n 



dr n dr \ d 



T 



(50) 



As the partition sums Z n , n > 1, cannot be expressed in terms of $3, or other 
special functions of mathematical physics, the internal energy, the heat capacity 
etc. cannot be written in compact form. However, due to the fact that Z n , 
n > 1, have similar behaviors near r = 0, (compare Eqs. (|48|) . (|49|) . we still 
shall have: 

c„~^^, (51) 

T 

which differs from the power low T-dependence which occurs for macroscopic 
and even microscopic systems. 

The temperature dependence of the heat capacity, for systems with n = 1, 
2, 3 and 4, is represented in Figure [TJ The heat capacity per particle is plotted 
with solid, dashed, dotted, dot-dashed line, for n = 1, 2, 3, 4 respectively. The 
maximum value of c n for n = 1, 2, 3, 4 occur at the following values of r: 

2.84, 1.16, 4.25, 5.2 . (52) 
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Their relative heights, defined as: 

Cn.max [T — 10) 

c„ (r = 10) 



(53) 



are: 

0.0376 , 0.2671 , 0.0321 , 0.0157 . (54) 

So, the effect manifests strongly at n = 2, but fades out for larger values of n. For 
large temperature, they show the same tendency as mesoscopic and macroscopic 
systems, i.e. tends to a horizontal asymptote, satisfying Dulong-Petit law. 

Two comments concerning the plots in Figure Q] should be done. First, it 
is quite unexpected that the maximum of the heat capacity is more visible for 
n = 2 than for n = 1. Second, we must be aware of the fact that Z n (max = 10) 
are good approximations for small t, so clearly more credible for t < 1 than for 
r~10. 

Another generalization of our initial one-fermion problem can be obtained 
considering, instead of a particle confined to a segment, the same particle con- 
fined to an elongated rectangle with edges a , b ; a <C b. The quantization along 
the axes Ox , Oy gives: 



k x — ~ i ky — — , (55) 



with the ground state one-particle energy levels: 

2 vn a 2 m o z 

The partition sum factorizes: 

Z = Z (q sx ) Z (q sy ) , q x = exp (—fiEi sx ) , q v — exp (—J3Eu y ) . (57) 

or: 

Z = \ [tf 3 (q sx ) ~ 1] [03 {Qsy) - 1] ■ (58) 
The heat capacity is a sum of contributions from each direction: 
C(q 

sxi Qsy ) = C (q sx ) + C (q sy ) , (59) 

where C (q sx ) is given by Eq (|39| . If we keep fixed one of the long edges of the 
rectangle, and curb the rectangle, superposing the other long edge on the fixed 
one, we get an empty cylinder; if we bend the cylinder, superposing its ends, we 
get an empty torus. For a particle moving on its surface, the quantization of its 
wavevector is similar to Eq. ([7]), and its thermodynamics can be obtained from 
the previous relations, with the index replacement s — > c. As q c — , we have 
again a situation when the heat capacity depends on the topology of the body. 
In principle, such a situation can be encountered in the case of a rectangle of 
graphene, and of a toroidal nanotube. The translational electronic contribution 
to the heat capacity for these two bodies will be different. If the heat capacity 
of such bodies can be measured with a precision of about 10~ 4 , this effect can 
be observed experimentally. 
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3 Heat capacity of a few-fermion system con- 
fined to various cavities 



Let us firstly consider a fermion in a cavity having the form of a rectangular 
prism, with square basis of area a x a and height b. We shall assume that the 
volume V — a 2 b is constant, and the cavity is described by the parameters a, 
b — V/a 2 . The movement of a particle inside the cavity is quantized according 
to the relations: 



7T , 7T ira 2 



k y = -, k z = - = — . (60) 



Putting: 



a 3 



(61) 



V 

the energy levels corresponding to each axis are: 

E\ x = Ely = , E\ z = E\ x e 2 , (62) 

2ma z 

E nx = E ny — Eis n 2 , E nz = e 2 E\ s n 2 . (63) 

If e <?C 1, the levels corresponding to the Oz axis are very dense, so they are 
give the dominant contribution to the heat capacity, but if e 3> 1 , they are very 
rare, and their contribution is negligible. The partition sum factorizes: 

ry rj, ^2 ry, S f Qx = (~ /3 El x ) = Q y 

Z = Z {q x ) Z (q z ) ,1 " (64) 

[q z = exp {-ps E lz ) = (q x ) 

or: 

Z=\ Wis (?«) - l] 2 [0a ((q.f) - 1] • (65) 

Due to this property, the heat capacity is a sum of contributions from each 
direction: 



C(q,s) = -k B p 2 



d 2 \nZ{q x ) , d 2 \nZ{q%) 



dp 2 d/3 2 
Let us put, similar to previous cases: 

1 E lx 1 E\ z e 2 



2C{q x ) + C{ql). (66) 



(67) 



t x k B T t z k B T k B T 

In the low-T limit, taking into account Eq. (I29p . we have the following limiting 
cases: 

Ife<l, 

C{q s s ) = C z ^C x . (68) 



1 




1 




< 




Tz 




T x 


1 




1 




> 








T x 



Ife>l, 

C(gS) = q,«C7 1B . (69) 
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In the high-T limit, Eqs. (® and {33J give: 



C (q, e) = | fc B + O (r 2 exp (- pi 2 r) ) , (70) 

so Dulong-Petit law is obtained, in its 3D form. 

If the prism is very long, a -C b or e <C 1, it can be bent until its opposite faces 
superpose each other, forming an empty ring; the quantization corresponding 
to the Oz axis is dominant, and the behavior of the system tends to that of a 
particle on a circle. 

The case of a particle in a spherical, cylindrical or toroidal cavity can be 
treated similarly to the ID problems discussed in Section 2. The situation 
is of course more complicated, for two reasons. First, the energy eigenvalues 
can be obtained only approximately, from the asymptotic values of zeros of 
Bessel functions (of half-integer order, for a sphere, and of integer order, for 
cylinder and torus). Second, the summation involves three quantum numbers; 
the azimutal one gives a trivial contribution, but the other two run over all 
integers, from 1 to oo. The final result is that the partition sum is not expressed 
any more in terms of $3, but in terms of more general functions: 

(q; a, (3)= J2 n a q {n+ ^ 

with a = 0, 1, 2 and (3 = 0, 1/2. We omit the detailed presentation of these 
functions, which are quite simple, but cumbersome and of small physical inter- 
est. The heat capacity of a particle confined to a sphere or a cylinder is quite 
similar to that presented in Figure [T] 

4 Conclusions 

The starting point of the investigations presented in this paper is the fact that, 
for few-particle systems confined to a small volume, their thermal behavior must 
be evaluated using a partition sum calculated discretely. In order to be able to 
treat the problem exactly, or at least in controllable approximations, only the 
free particles (in fact fermions) without internal structure have been considered. 
In this case, only the translational degrees of freedom contribute to the partition 
function. According to the remark just made, the partition function of a small- 
volume system should not be evaluated as an integral over momenta, as usually 
done for large systems, but as a sum over the quantized values of the particles 
wave vectors. 

As recently discussed by Toutounji [1], if we evaluate in this way the canoni- 
cal partition sum for an one-fermion system, confined to a segment, and compute 
the heat capacity using the standard thermodynamic formulas, we find that this 
quantity shows a maximum at low temperature. Even if, with increasing T, the 
heat capacity tends to a horizontal asymptote, its low-T behavior differs quali- 
tatively of that of a macroscopic body, which increases smoothly from zero, at 
T = 0, to a constant value, at large temperatures, according to the Dulong 
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- Petit law. In fact, the same macroscopic behavior occurs, at least qualita- 
tively, for mesoscopic systems too, as shown recently by Lungu, who considered 
systems of 14, 76 and 820 particles [2J. 

In this respect, our paper reveal two new aspects. First, it shows that the 
heat capacity falls at T = according to an exponential law, merely than a 
power law, as in mesocsopic and macroscopic systems. Second, the maximum 
of the heat capacity, discussed by Toutounji for a system with n = 1 particles, 
is even stronger for n = 2, but fades out for n — 3 and n = 4. 

Using the same theoretical methods, we calculated the heat capacity of few- 
body systems having the same local geometric properties, but different topology 

- for instance, a segment and a circle with identical length, or a cylinder and a 
torus obtained by bending this cylinder - and we found that the heat capacity 
depends on the topology of the body. Also, few-particle systems confined to 
cavities having the same volume, but different geometries are different. Even if 
these conclusions might be surprising at first sight, they are simple effects of the 
sensitivity of the thermal properties of few-particle systems to the specificities of 
energy spectrum, at low T, and do not contradict the well known results valid 
for macroscopic bodies, obtained theoretically after the thermodynamic limit 
(N — > oo, V — > oo) is taken. 

In the same time, these results are simple and pedagogical illustrations for 
the peculiar behavior of nanoscopic systems, as compared to macroscopic, and 
even with mesoscopic ones. These effects - for instance topology-dependent heat 
capacity - could be observed experimentally, in systems composed of benzene 
molecules or nanotubes. 
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Figure 1: The heat capacity per particle for a system of n fermions; n = 1 : 
solid line; n = 2 : dashed; n = 3 : dotted; n = 4 : dot - dashed. 
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